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ABSTRACT 

The Cold Dark Matter (CDM) model for galaxy formation predicts that a significant 
fraction of mass in the dark matter haloes that surround L ~ L* galaxies is bound 
in substructures of mass 1O 4 -1O 7 M0. The number of observable baryonic substruc- 
tures (such as dwarf galaxies and globular clusters) falls short of these predictions 
by at least an order of magnitude. We present a method for searching for substruc- 
ture in the haloes of gravitational lenses that produce multiple images of QSOs, such 
as 4-image Einstein Cross lenses. Current methods based on broadband flux ratios 
\q • cannot cleanly distinguish between substructure, differential extinction, microlensing 

C^) and, most importantly, ambiguities in the host lens model. These difficulties may be 

£N| , overcome by utilizing the prediction that when substructure is present, the magni- 

' fication will be a function of source size. QSO broad line and narrow line emission 

regions are approximately ~ lpc and > 100 pc in size, respectively. The radio emis- 
sion region is typically intermediate to these and the continuum emission region is 
I 1 much smaller. When narrow line region (NLR) features are used as a normalisation, 

the relative intensity and equivalent width of broad line region (BLR) features will 
respectively reflect substructure-lensing and microlensing effects. Spectroscopic ob- 
servations of just a few image pairs would probably be able to cleanly extract the 
desired substructure signature and distinguish it from microlensing - depending on 
the actual level of projected mass in substructure. In the rest-optical, the H/3/[0iii] 
k> region is ideal, since the narrow wavelength range also largely eliminates differential 

reddening problems. In the rest-ultraviolet, the region longward and including Lya 
' may also work. Simulations of Q2237+030 are done as an example to determine the 

level of substructure that is detectable in this way, and possible systematic difficulties 
are discussed. This is an ideal experiment to be carried out with near-infrared integral 
field unit spectrographs on 8-m class telescopes, and will provide a fundamentally new 
probe of the internal structure of dark matter haloes. 
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1 INTRODUCTION 

Galaxy formation models cast in a ACDM universe have ex- 
perienced considerable success to date. With a fixed choice 
of cosmological parameters, it is possible to make specific 
predictions, for example based on N-body simulations, re- 
garding the structure of individual dark matter haloes. One 
of the challenges in this field is to calculate the level of the 
continuously-infalling substructure that will survive within 
large haloes. Present N-body simulations indicate that there 
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should be plentiful substructure in a L, galaxy's halo, with 
masses above 1O 7 M0; below this scale limitations in reso- 
lution becomes important. If these subhaloes contain stars, 
there should be signi ficantly more dwarf galaxies near th e 



Galaxy than are seen (Moore et al. 1999; Klypin et al. 2001) 



Various solutions to this excess problem have been proposed, 
including warm d ark matter (which erases substructure on 
small scales; e.g. Avila- Reese et al. 2001); self- interacting 



dark matter (which causes subhaloes to evaporate; Spergel 
& Steinhardt 2000); and radiation-driven ionization of the 
intergalactic medium such that star formation in dwarf- 
size galaxies is su p pressed at early epochs, as discussed i n 
Somerville (200$ , [Bullock, Kravtsov fc Weinberg (200l[ ), 
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and |Klypin et al. (1999| ). If the substructure does exist, its 
gravitational influence should be detectable, and therefore 
this prediction of ACDM in particular may be tested. There 
is already some evidence in support of the need for substruc- 
ture to explain the relative magnifications in strong lenses; 



disenta lgling its signature from other (equally interesting) 



effects is the scope of this paper. 

From the positions of the images and lens galaxy of 
a 4-image lens system, a parametric smooth mass model 
for the lensing galaxy and halo can be computed (including 
external shear), along with predictions for the three inde- 
pendent flux ratios. Measuring deviations between the pre- 
dictions and the (differential-reddening-corrected) observed 
integrated fluxes in real Einstein Cross lenses, provi des a di- 
rect est imate of t he surface d e nsity in substructure ( Mctcali 



& Mad ku 200 j; |Chiba 2002| ; [Metcalf fc Zhao 2002| ; palal 
& Kocl anek 2002). Substructure has long been considered 



a possible explanation f or inconsistencies in th e observed 
flux ratio of B1422+231 (|Mao fc Schneider 199^). However, 



this approach has a vital flaw. The predicted flux ratios can 
be strongly dependent on the parametric lens models used, 
making a ny measurement of the ma ss fraction in substruc- 



ture (e.g. Dalai & Kochanek 2002) suspect. The observa- 



tional technique proposed here largely avoids this impor- 
tant problem, as well as the problems of intrinsic variability 
(on the timescale of the typical time-delay between images) , 
and the effects of microlensing, with the possible change of 
continuum slope as well as intensity. 

In §P the main features of substructure and its lens- 
ing effects are summarized. Simulations demonstrating the 
i mpo rtance of source-size on magnification are presented in 



2.3 



for a specific lens, Q2237+030. The spectroscopic ap- 
proach proposed here is discussed in with the most dif- 
ficult problems and their possible solutions addressed. In 
§^| we discuss selection criteria for lens candidates for this 
experiment, as well as modes of observations that would 
be sufficient for obtaining results. The main conclusions 
are summarized in In the Appendix an estimate of 
the observational uncertainties is calculated. All cosmo- 
logical calculations are carried out for a ACDM cosmol- 
ogy, with = 0.7, Qo = 0.3, though none of the argu- 
ments (except the level of substructure one might expect) 
depend strongly on this choice. The Hubble parameter is 
Ho = 65/i 6 5 km s" 1 Mpc~\ 



2 SUBSTRUCTURE & GRAVITATIONAL 
LENS MODELS 

2.1 substructure and small-mass dark matter 
haloes 

The highest-resolution N-body simulations of structure for- 
mation are curren tly limited to ~ 1O 6 M0 particles (e.g. 
Moore et al. 1999). Down to the scale of 10 particles (~ 



10 Mq) there is clear evidence for surviving substructure 
within large dark matter haloes with a power-law mass func- 
tion, dN/dm oc m ~ a , where a ~ 1.9 1. For observed dwarf 
galaxies, a ~ 2.35 ( Klypin et al. 2001 ). In these simulations 



there is a continuous process of infall and merging work- 
ing against tidal stripping and dynamical friction to main- 
tain the substructure population. Within the virial radius 



of a galactic halo, about 10-15% of the mass is in structure 
of mass Si 1O 7 M0. Analytic arguments suggest that future 
simulations with a larger dynamic range will find that sub- 
structure survives down to far smaller masses than are cur- 



rently accessible by direct computation (Metcalf & Madau 



2001) 



In addition to the substructure inside the primary lens, 
there will also be small haloes in intergalactic space in front 
of and behind the primary lens. Zhao & Metcalf (2002; in 
prep.) show that the importance of these two populations 
to compound lensing is approximately equal, if the tidal de- 
struction of haloes within the lens is negligible. Tidal de- 
struction will increase the relative importance of intergalac- 
tic haloes although the importance of this cannot yet be 
accurately quantified. Despite this, in our estimates and sim- 
ulations we treat all the haloes as if they are at the same 
redshift as the primary lens. Since the strength of a lens of 
fixed mass reaches a maximum when the lens is half way 
to the source, measured in angular size distance, and the 
primary lens is usually near this point, an estimate of the 
surface density in substructure under this simplification will 
be an estimate of the minimum density required to produce 
the observed lensing effects. 

The importance of haloes to lensing is strongly depen- 
dent on their internal structure. To quantify the internal 
structure in some of the estimates given below, we use the 
parameter 7(7-) = din M(r) /dlnr, where M(r) is the mass 
of the halo within a projected radius of r. The tidal radius of 
a subhalo is set by the requirement that the average density 
within that radius be proportional to the average density of 
the host halo within the orbit of the subclump - the con- 
stant of proportionality depends on the structure of the host. 
This means that the average density p c = 3m/47rr(m) 3 , 
where r(m) is the radius of a subhalo, will be independent of 
mass. Additionally, in the hierarchical structure formation 
model, the initial density of a halo is proportional to the 
average density of the universe at the time when the halo 
first collapsed. In the CDM model all the small-mass haloes 
of interest here are almost coeval. This results in the same 
conclusion for the intergalactic haloes as for the subhaloes - 
p c is independent of mass although it may have a different 
value for the two populations. We will take p c to be a free 
parameter in our discussions. 

2.2 the influence of source size on lensing 

We can make some estimates for how compound lensing 
will be affected by the size of the source QSO using the 
simple model of substructure outlined above. Let us brake 
the magnification matrix, Ay, into a contribution from the 
smooth lens and a perturbation caused by the small clumps 
us say a clump will have a significant 
influence on the magnification if <5Ay is of order e or larger. 
The value of e will depend on the magnification caused by 
the unperturbed smooth lens (Sfi/fj, ~ [i e where \jl is the 
smooth lens magnification). The magnification of a source 
with finite size is the point-source magnification integrated 
over the area of the image. Because of this, a finite source 
sees a magnification that is smoothed on the scale of its 
own size. For a halo to have an influence, the region where 
its contribution to the magnification matrix is greater than 
e must be larger than the image of the source. If we take 
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Figure 1. The minimum mass of a small halo that can have an 
influence on the compound lensing as a function of the source 
size. The solid line is from Eqn. ^. The dashed line is Eqn. ^ with 
R = 1 kpc and cr h = 200 km s — 1 . In both cases hes = 1 and 
e = 0.1. If the magnification is high (fi X, 10), e = 0.01 would be 
more appropriate. The general size range for different emission 
regions of a QSO are labelled (BLR - broad line region; NLR - 
narrow line region). The central accretion disc of a QSO which is 
believed to be the source of the continuum emission has a scale 
of 10 -4 - 10 -3 pc; off the scale to the left. 



the physical size scale of the source to be £ s , this require- 
ment puts a lower limit on the mass of haloes that can be 
important for lensing: 



£7rE c 

1 + 7 



/3 -| 3/(3-7) 



(1) 



where the critical density is E c = (4ttGDiDi s /c 2 D s ) _1 and 
the angular size distances to the lens, source, and from the 
lens to the source are denoted by Di, D s , and D\ B , respec- 
tively. The logarithmic slope 7 is evaluated at the radius 
where 8A(r) — e. 

It is useful to put some typical numbers into Eqn. |l| 
Consider a source at z s — 2 and a lens at zi = 0.5. In this 
case E c = 1.93 x 10 9 /!,6 5 M Q /kpc 2 and Di/D s = 0.73. For a 
point-mass (the most condensed possible halo), 7 = and 



m c ~ 3.2 x 10 3 e/i 65 ' 

' PC 



Ma 



(point-mass lens). (2) 



For less condensed haloes we must make an assumption 
about their internal density, p c . The characteristic density of 
the intergalactic haloes, p c , is set by the process of structure 
formation at high redshift. If these haloes fall into a host lens 
they will remain intact until the average density of the halo 
within their orbit, ~p h (R), becomes of order the density of the 
subhalo, then it will be tidally stripped. If all the subhaloes 
are tidally truncated, their density can be calculated from 
the average interior density profile of the host, ~p h (R). If 
the host is modeled as a singular isothermal sphere (SIS), 
p c = 2p h (R). This results in a cutoff mass of 



, . i , , / , £ a \ 200 km R ,„ T _> 

m c ~1.7xl0 4 M Q e/i 65 — = (SIS), 

pc I a h s kpc 



where is the velocity dispersion of the host lens. The 
haloes could be denser than assumed here if they were born 
denser, or they could be less dense if they were born less 
dense and primarily reside outside of the primary lens. Since 
no accurate simulations of structure formation at these small 
mass scales have been done, we will take the above estimates 
as benchmarks for our discussions. 

The estimates of m c are plotted in Figure |l| as a function 
of source size along with the approximate size scales for dif- 



ferent emission regions of a QSO as discussed in section 3.1 



The image of a source of a given size will be affected by halo 
masses above this limit. For example, the image of a 1 pc 
source will reflect the influence of haloes with m Si 1O 3 M0 
while the image of a 100 pc source will only be affected by 
haloes with m Si 1O 6 M0. For source sizes below ~ 10 -2 pc 
(such as the accretion disc of a QSO), m c is less than a solar 
mass, making microlensing by ordinary stars possible. Fig- 
ure H shows m c as a function of lens redshift. The drop in 
m c near zi — in Figure |^ for the point masses is somewhat 
misleading here because the Einstein radius becomes very 
small near z\ = making these cases extremely rare. The 
value of m c for any realistic substructures lies somewhere 
between its value for point mass lenses and SIS lenses. 

This mass cutoff alone indicates that the magnification 
should be a strong function of source size if there is a large 
population of haloes with m ,$ 1O 7 M0 . There are other fac- 
tors that influence how magnification will depend on source 
size. Larger sources are more likely to be eclipsed by rarer, 
larger haloes. On the other hand, if the haloes have high 
density cusps in their centers - as expected from simula- 
tions - smaller sources can reach higher magnifications. To 
determine accurately how the magnification distribution will 
depend on source size and halo mass distribution requires 
numerical simulations. 



2.3 lensing simulations: Q2237+030 

We take the image positions and lens galaxy position of the 
4-image QSO lens system Q2237+030 and fit a elliptical 
power- law plus shear model to them for the smooth lens 
component (see 



Metcalf fc Zhao 2002 for details). To cal- 



culate the magnifications with substruc ture, an improved 
versio n of the ray tracing code used by 



Metcalf & Madau 



(2001) is used. The lensing equation is solved to obtain the 



(3) 



source position for each point on a coarse grid surrounding 
the unperturbed image position. The closest grid point to 
image position is found and a new smaller grid is set up 
around that point. This process is repeated until the image 
takes up a large fraction of the grid and the desired accu- 
racy is reached (usually Sfi/fi — 10 -3 ). This grid-refinement 
technique allows us to achieve the large dynamical range re- 
quired to calculate the magnifications of sources with sizes 
ranging from less than 1 pc to several hundred pc while keep- 
ing the substructure fixed. 

The substructure is modeled as singular isothermal 
spheres truncated such that ~p c — 6.9Mq pc -3 and with a 
mass function of dN/dm cx m -2 . The total surface den- 
sity in substructures with 10 3 Mq < m < 10 6 M© is fixed 
at 10% of the total surface density according to the smooth 
model. Random realizations of the masses and positions of 
the subclumps are made and the magnifications with several 
different sources are calculated in each one. 
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Figure 2. The minimum lensing mass as a function of the source 
size as defined in Eqn |lj. The solid curves are for point mass 
substructures and the dashed are for singular isothermal spheres. 
The lower of each set of curves are for a source size of l s = 1 pc 
and the upper for i s = 1 pc. The curves that stop at 2; = 1 are 
for source redshift z s = 1. The others are for z B = 2. The SIS are 
tidally truncated with R = 1 kpc. A realistic case would probably 
lie somewhere between the point mass and the SIS curves of each 
type. In all cases, h^5 = 1 and e = 0.1. 



Figure 4. The differences in magnification ratios between image 
A and image D of quad lens Q2237+030 using three different 
source sizes (lpc, 10 pc and 100 pc). A smooth, substructure-free 
lens would predict that the three magnification ratios would be 
very close to equal, the center of the plot. These are the same 
1 000 realizations as in Figure jj| One quarter of the realizations 
have |Ami4(10 pc) — Ami4(l pc)| > 0.05 mag and one quarter 
of them have |Ami4(l pc) — Ami4(100 pc)| > 0.18 mag. 
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Figure 3. Simulated magnification ratios between image B and 
image D of quad-lens Q2237+030 using two different source sizes. 
Each circle is a different realization of the substructure. There are 
1 000 realizations shown. Without substructure, A7TI24 should not 
depend on source-size in this range so all the points would be on 
the 45° line shown. The range of ratios along the 45° line is larger, 
but excursions in this direction can be degenerate with the smooth 
lens model. Excursions perpendicular to the line unambiguously 
indicate small-scale structure. 



Figure [| shows the flux ratio, measured in magnitudes, 
of image A to image D (marked '1' and '2') in Q2237+030 for 
1 000 realizations of the substructure. The ratios for source 
sizes of lpc and 10 pc are shown. If there were not a sub- 
stantial number of haloes with mass less than ~ 10 6 Mq, we 
would expect these two ratios to be equal and the circles in 
the plot to he on the diagonal line. No smooth model can 
account for points being off this line. 

Figures H through ]a show simulations of all three mag- 
nification ratios for Q2237+030. In each case, three differ- 
ent source sizes are used. If the lens were smooth all three 
magnification ratios would be the same and all the dots in 
the figures would be located precisely at the origins in these 
plots. Figure shows the cumulative distribution of the sim- 
ulated differential magnification ratios. We can see that in 
half of the cases the differential magnification ratio is larger 
than 0.1 mag. If this level of measurement accuracy can be 
attained, only 1-3 lens systems will be required to rule out 
this level of substructure. 

Our model for Q2237+030 predicts magnifications that 
are relatively modest (;xa = 2.28, fiB = —0.939, [ic = 2.047, 
[id = —2.005). Since the substructure has an influence on 
the magnification that is proportional to the smooth compo- 
nents magnification (Sfj, oc we expect that these results 
are fairly conservative. Some published models for other 
quad lenses predict magnifications as large as 100. 

We can see from these simulations that it is possible to 
detect substructure, or intergalactic dark haloes, in a way 
that is independent of any lens model, if we can measure 
the flux ratios to ~ 0.1 mag for sources that differ in size 
by a factor of ten and are centered on the same object. 
This can be done with a very small number of lens systems. 
Interpreting the results in terms of the fraction of mass in 
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Figure 5. Same as in Figure |4|, but for images B and image D 
of Q2237+030. In this case one quarter of the realizations have 
|Am24(10 pc) — Am,24(l pc)| > 0.05 mag and one quarter of them 
have |Am 2 4(l pc) - Am 2 4(100 pc)| > 0.07 mag. 



Figure 7. Here is a histogram the cumulative distribution of 
the differential magnification for the three independent ratios in 
Q2237+030. On the y-axis is the fraction of simulations out of 
1 000. The source sizes used are 1 pc and 100 pc to reflect the 
approximate sizes of the BLR and the NLR, respectively. 
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cantly over the time span of several years (Kaspi et al. 1996) 



Figure 6. Same as in Figure ^ but for images C and image D 
of Q2237+030. In this case one quarter of the realizations have 
|Ar?i34(10 pc) — Am,34(l pc)| > 0.05 mag and one quarter of them 
have |Am,34(l pc) — Arrt34(100 pc)| > 0.17 mag. 



small mass haloes and the slope of the mass function will 
require some modeling of the smooth lens component, but 
even then the results will be much less model dependent 
than the more traditional method based on single broad- 
band measurements of the flux ratios. 



3 THE SPECTROSCOPIC APPROACH 

3.1 the internal structure of QSOs — a sketch 

In QSO spectra, the 'narrow line region' (NLR) emission 
lines (e.g. [O in] AA4959, 5007) do not vary in flux signifi- 



For this reason, they have been used for calibration in long- 
term spectroscopic monitoring observations of low redshift 
QSOs for reverberation mapping experiments, by which the 
time delays between features of different ionization states 
may in principl e be used to m ap all the components that 
make up QSOs ( Peterson 1993 ). Photo-ionization arguments 
and cases where the NLR-emitting region is imaged directly, 
indicate that the size of the NLR extends out to 100-1000 pc. 

In the full unification picture for active galactic nuclei 
(AGN) including QSOs, the variability timescales of dif- 
ferent spectroscopic features, and the time delays between 
them, are used to map out the size and structure of dif- 
ferent componen ts with very different characteristics (c.f. 
Antonucci 1993). The basic idea consists of a supermas- 
sive black hole (Mth > 1O 6 M0) in the core, surrounded by 
an accretion disc, which produces the nearly flat-spectrum 
continuum light seen. The continuum flux varies on very 
short timescales, less than a day, and so must be very com- 
pact (around 100 AU or ~ 5 x 10 _4 pc). It is also known 
to be microlensed in the case of Q3327+ 030 which puts 
an upper limit of 2 000AU on its size (Yonehara 2001, 
and references therein). Beyond that, there are permitted 
lines such as the Balmer-series lines that are very broad 
(dfwhm > 10 3 kms _1 ) due to gravitationally- induced mo- 
tions. The characteristic size of this 'broad line region' 
(BLR) apparently scales with the intrinsic luminosity of the 
host QSO (and therefore with the m ass of the central black 
hole; Kaspi et al. 1996| Wandel 1999), such that in luminous 
QSOs (as opposed to, e.g., Seyfert l's), the size of the BLR 
is on the order of six light-months (or ~ 0.3 pc). According 
to long-term studies of QSOs, the typical (rest-frame) flux 
variations in the QSO continuum are on the order of 10- 
70% (tho ugh occasional flu ctuations as high as ~ 50 x are 
possible; Ulrich et al. 19971), while the BLR variations are 



smaller by a factor of 2-4 ( Maoz et al. 1994 ) . In the observed 
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Figure 8. In the rest-optical, the shown rest-waveleng th region 
(from the SDSS QSO template, Vanden Berk et al. 2001 ), includes 



the NLR [Om] AA4959, 5007 and the BLR H/3 emission features. 
Where the NLR emission line is used for normalising the spectra 
of each lensed image, any residual in the relative intensity of the 
BLR H/3 emission lines after accounting for the relative magnifi- 
cation in the relevant gravitational lens model, will most likely be 
due to compound-lensing due to substructure in the foreground 
'smooth' lens gravitational potential, whereas intensity variations 
in the continuum will be due to a combination of substructurc- 
and micro-lensing, and other effects as discussed in the text. 



frame, these timescales are stretched by an additional factor 

Of (1+Zs). 

There is some evidence for anti-correlations between the 
rest- ultraviolet continu um flux and th e equivalent width o f 
certain emission lines (Baldwin 1977, Green et al. 2001), 
which may also be related to the correlation bet ween the 
cent ral black hole mass and the continuum flux (Wandel 
general trend claimed, is that the more lumi- 
nous QSOs have less pronounced broad-line emission line 
components, the [Om]/H/3 central-flux ratio is small, and 
the rest-ultraviolet and -optical Fe n multiplet emission lines 
are more prominent, obfuscating the 'true' local continuum 
level that underlies these emission lines. There is strong 
(brightening) evolution with redshift observed in the lumi- 
nosity function of QSOs, so that in general, high-redshift 
optically-selected QSOs are less likely to have strong [O in] 
emission. The QSO templates from the Large Bright QSO 
Survey (LBQS; Francis et al. 1992) and the Sloan Digital 
Sky Survey (SPSS; |Vandcn Berk et al. 2001| ) are compiled 
from objects over a broad range of redshifts, relatively few at 
z < 1 (where the rest-optical window around H/3 is still vis- 
ible). Therefore, the H/3 region in these templates is drawn 
from objects that are relatively local and of lower luminosity, 
and so generally show a misleadingly large [O in] /H/3 ratio, 
if they are to be applied to higher-redshift considerations. 

The only definite way to determine the characteristics 
of a particular system is by direct observations in the rest- 
optic al. For example, in th e case of B1422+231 at z s = 
3.62, |Murayama et al. (1999| ) find that the Fe n AA4434-4684 
emission is relatively weak compared to H/3, similar to lower- 
redshift QSOs. 



3.2 the spectroscopic approach for substructure 

Consider a case where the spectrum of all images in an Ein- 
stein Cross are obtained, in a wavelength range where both 
NLR and BLR emission features are visible at high signal- 
to- noise (SNR» 10), and with sufficient spectral resolution 
to clearly resolve the shapes of the lines (v res ~ 200 kms^ 1 ). 
The NLR emission lines are used to normalise the spectra 
to each other, and we wish to attempt to disentangle the 
substructure signature from other effects. There are several 
systematic errors to consider. 

Differential Reddening — As each lensed image fol- 
lows a different line of sight through the lensing galaxy, if 
the gas and dust in the galaxy are sufficiently patchy or 
position dependent, it is possible that there will be dif- 
ferential reddening that may affect the meas ured contin- 



uum slope and relative emission line fluxes. Falco et al 
(1999} ) measured the A E(B - V) in 23 lens galaxies us- 
ing HST broadband colours, under the (optimistic) assump- 
tion that variability and microlensing were not important. In 
recent broadband-magnitude-based substructure estimates, 
the colour corrections estimated by the above work were 
applied, which, in light of the possibility of microlensing, 
may introduce unquantifiably systematic errors for any one 
system. In the spectroscopic approach, choosing NLR and 
BLR emission lines that are close in wavelength, the differ- 
ential reddening is very small. Taking the media n value of 



A E(B -V)= 0.04 mag in the |Falco et al. (1999| ) optically- 
s elected sample, and assuming a MWG-type extinction law 
( |CardeUi et al. 198S| ), we estimate that flux ratio between 
[O in] AA4959, 5007 and H/3 for z 3 = 3 and zi = 1 will vary 
by less than 0.07% and can therefore be ignored. For com- 
parison, between [O n] A3727 and Mg n A2800 the variation 
will be about 1.4%. 

Microlensing — Some microlensing by stars proba- 
bly occurs in m ost lenses, but it is most cl early seen in 
quadruple lenses (Witt, Mao & Schechter 1995, and citations 



thereof). Microlensing events last months or years, much 
longer than the time delay between images, but (especially 
at caustic crossings, which can occur quickly) c an magnify 



an image by as mu ch as a whole magnitude (e.g. Witt, Mao 
&i Schechter 1995). For super-microarcsec sources, like the 
BLR of a QSO, stellar microlensing is insignificant (Fig- 
ure [l]), but for smaller sources, like the more compact QSO 
non-stellar continuum, microlensing will affect the magnifi- 
cation. Hence, observations of the relative flux of the broad 
component of BLR lines in the images (e.g. after any narrow 
component has been fit out) will reveal the effects of lensing 
by CDM substructure alone, while their equivalent widths 
can be used to probe ongoing microle nsing because they ar e 
measured relative to the continuum (McLeod et al. 1998). 
Thus, in a spectroscopic experiment as we have described, 
it would be possible to place a limit on microlensing and 
substructure lensing simultaneously. 

If the QSO accretion disc has a steep temperature gradi- 
ent, it is also possible for microlensing of the disc to produce 
wavelength-dependent magnification in the co ntinuum of the 



affec ted image (as seen in HE 1104-1805 by Courbin et al 
200C). This does not affect the BLR emission line fluxes, as 
compared against NLR lines, though it is also spectacular 
evidence for microlensing. 
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Intrinsic QSO Variability — The continuum flux 
variations in QSOs can be quite large and rapid, on 
timescales comparable to the time delay between images. 
The BLR emission line fluxes vary less dramatically, and 
only reflect the more pronounced fluctuations seen in the 
continuum, with a several-month time-lag in the rest frame. 
For source redshifts z 3 m 2, the observed-frame time-lag 
becomes close to a year. As long as variations in the ob- 
served frame happen on timescales that are significantly 
longer than the time delay between images in a given lens 
system, using the BLR lines as a substructure probe should 
be robust. 



4 DISCUSSION 

The ideal target is a four-image lensed QSO with an excep- 
tionally well-determined (smooth) gravitational lens model, 
at a redshift such that the ideal H/3/[0 in] features are red- 
shifted into convenient observable windows, unobscured by 
strong atmospheric emission bands. It is also desirable for 
the target QSO to have a high [Om]/H/3 ratio, which typi- 
cally also means that the H/3 emission will be dominated by 
its BLR component, and also that the Fen complex short- 
ward of H/3 will be weak, such that the underlying local 
continuum can be measured securely. The redshifts of the 
sources in most Einstein Crosses are known; their H/3/[Oiii] 
features have only been measured, at low signal-to-noise, in 
very few cases. Therefore, in many cases exploratory obser- 
vations may be necessary. Several systems will need to be 
studied in detail in any case, to build up the statistics needed 
for a definitive measurement. 

As illustrated in Figure ^, in the rest-optical the 
[O in] AA4959, 5007 and H/3 lines are close, but separated 
enough to be distinguishable at even moderate spectral res- 
olution. At the redshifts of most gravitationally lensed QSOs 
(typically z > 2), these features are shifted into one of the 
near-infrared windows at A > l/im (JHK). 

In the Appendix we investigate the expected statistical 
uncertainties in a measurement the ratios between lines and 
contin uum. Based on a 'typica l' QSO (e.g the SDSS tem- 
plate; Vanden Berk et al. 2001), the BLR to NLR emission- 
line ratio may be measured with an error of ~ 0.3/ SNR, 
where SNR signal to noise ratio per resolution element 
with which the continuum is measured. Likewise, the con- 
tinuum to NLR emission-line ratio can be measured with 
error a c ~ 10/ SNR. If in the continuum a SNR ~ 100 is 
achieved, then based on our simulations in {2.3 (with ~ 10% 



of the surface density in compact objects below ~ 1O 6 M0, 
as seems to be the case in ACDM), then we expect the line 
ratios between the images of a typical lens system to differ 
by an amount larger than the measurement error. The situa- 
tion is actually a bit better than this because all the images 
can be used to fit the full line shapes simultaneously. Ob- 
servations of additional lens system can be used to improve 
the sensitivity to the substructure mass fraction. 

These levels of the SNR are easily feasible with 8-m class 
telescopes for many Einstein Crosses and relatively moder- 
ate amounts of observing time. For instance, the Einstein 
Cross Q2237+030 QSO images, at z 3 = 1.69, have infrared 



magnitudes of H » 16.0 (from CASTLES[|). With 3 hours 
of integration on an 8-m telescope with v re3 ~ 50kms _1 , 
the relevant part of the continuum, will be detected with 
SNR ~ 100. 



5 CONCLUSIONS 

In ACDM, a significant fraction of the dark matter in L ~ 
L* galaxy haloes is in coherent substructure, with masses 
of 10 4 -10 8 Mq , which does not seem to have baryonic, or 
at any rate luminous, counterparts. The substructure and 
small-scale intergalactic structure should be detectable by 
its gravitational signature in multiply-imaged QSOs. If such 
substructure is sufficiently abundant along the images' lines 
of sight, the light from different emission regions of a QSO 
will be magnified differently. This would be a model inde- 
pendent signature of small scale substructure. 

Broadband photometry is primarily a measurement of 
the continuum flux, which is affected by rapid fluctuations 
intrinsic to the QSO, and by microlensing caused by stars. 
The magnitudes of both of these effects are at least com- 
parably to what substructure is expected to produce, and 
conventionally can only be tracked by extensive monitoring. 

We propose that these fundamental limitations can be 
largely overcome by obtaining high signal to noise spectra 
(SNR ^ 100) of two or all lensed images simultaneously, 
that include prominent NLR and BLR features. The NLR 
features are used to normalise the spectra to each other, and 
then the relative fluxes of the BLR line(s) should be primar- 
ily due to substructure differential magnification. Intrinsic 
variability in BLR lines is observed to amount to ~ 10-35%, 
but relatively smoothly, and on characteristic timescales of 
greater than six months in the rest frame (and l+z 3 times 
longer in our observed frame). Compared to the time delay 
of only days or weeks between the images in typical Ein- 
stein Crosses, the intrinsic BLR fluxes can be assumed to 
be constant. 

The ideal setup is with an integral field unit (IFU) in 
the near- infrared (NIR), targeting the rest-optical [Oin] & 
H/3 emission lines. With the IFU setup, it is possible to get 
high-quality spectra of all images (and possibly of the lens, 
as well) , simultaneously. The alternative of a set of longslit 
observations can be made more efficient by aligning on pairs 
of lens images, positioning the slit close to the parallactic an- 
gle to minimize the effects of differential refraction. Efficient 
implementation will require an 8-m class telescope, with ex- 
cellent seeing conditions (~ 0.6"), perhaps with adaptive 
optics. 

This experiment will provide a concrete measurement, 
or at least place a severe upper limit on the fraction of dark- 
matter substructure in galaxy haloes and intergalactic space 
as predicted by the popular CDM models of structure forma- 
tion. These models are presently in a state of crisis because 
of the lack of observations of this substructure. 



The CfA- Arizona Space Telescope LEns Survey of gravitational 
lenses, at cfa-www.harvard.edu/castles/ 
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Parameter Values Estimated Errors in units of a a 

R+ AAt ui c acC(X)i A n A b cr„ a c a a /a cr A „/A n a Ab /A b 

3000 50,000 1.0 22.2 1.7 x lO" 16 300 5,000 0.255 13.1 0.538 0.531 1.19 
3000 100,000 1.0 19.9 1.7 x lO" 16 500 10,000 0.148 6.88 0.309 0.308 0.746 

* Resolution of spectrograph, R = A /AA. 
t A A, A n and A;, are expressed in km/s. 
tThe continuum flux in erg cm - 1 s — 1 A 1 . 

Table X. 



APPENDIX 

EXPECTED UNCERTAINTIES IN LINE STRENGTH RATIOS 

We wish to estimate how well the ratio of the strengths of the broad and narrow lines can be measured. To do this we will 
model the flux in the ith resolution of the spectrum as consisting of four contributions: the broad lines (cj'PVj (Aj,)), the 
narrow lines (/ n VK™(A n )), the continuum emission (c'd) and instrumental noise (Ni). The functions Wf, Wp and d are 
normalised so that their sum over all resolution elements is one and their widths (and perhaps some other shape parameters) 
are Aj, and A n . We are interested in the ratio uj = uj' / f n . The flux is then 

Ff = hriuWfiA^ + wriA^ + cCJ + Ni (Al) 

where /; is the extinction and c = c'// n . We will rewrite the extinction curve as i;/ n = ahi with the normalisation ^ . hiW™ = 
1 so that a is the total flux in the narrow lines after extinction correction. If we assume Gaussian, uncorrelated noise the log 
of the likelihood function will be 

- InC = J2 ^2 [o/n(wWf (A*) + W?(A n ) + cd) - Fi\ 2 + ln(2 7 r(T ! 2 )/2. (A2) 



For our present purposes it is convenient to replace the sum over resolution elements in (A2) with an integral over wavelength 
by making the substitution J^. — » S\~ J d\ where 5\ is the width in wavelength of the resolution elements. We will assume 
that the noise is constant over the wavelength range considered and given by a n . 
The Fisher information matrix is defined as 

92 (A3) 



dpidp2 

where pi and P2 are two parameters that are to be constrained by the measurements. The variance of the minimum variance 
unbiased estimator of a parameter p is given by the Cramer-Rao bound 

4 = (F-\ P . (A4) 



To calculate this we must identify the parameters in (A2) that are to be fit to the spectrum. The region of the spectrum 



we are interested in, shown in Figure M, is relatively small c ompared to the complete spectrum and it is in region where the 



continuum has a shallow minimum in template QSO spectra( Vanden Berk et al. 2001). For simplicity we will approximate the 



extinction, /i(A), and continuum, C(A), as constant over the region of interest, but we will allow their normalisations to vary. 
Formally, this is a valid approximation if d In C (A) /dX, d In h{ A) /d\ <C 1 / AA where AA is the range of wavelength considered. 
We will also assume that the redshift of the QSO is well established so that the positions of all lines are known. This leaves 
five parameters to be determined - u>, c, a, A n and A;,. There are then 15 independent elements in the Fisher matrix which 
we calculate. In general the profiles of the lines are functions of multiple parameters which could be added to the list, but 
for simplicity we will leave these out here and take the lines to be Gaussian. Once these matrix elements are calculated the 
matrix can be inverted numerically. 

It is convenient to express error estimates in terms of the signal to noise per resolution element with which the continuum 
can be measured. In our parameterization this is 

- v^AX 

(Jo = — jt-. (A5) 

A noise level of o ~ 0.01 may be reasonable. Table |l| shows the estimated errors for several choices of fiducial parameters. 
The ratio of the continuum flux within the range AA to the narrow line flux, c, is estimated from the height of the narrow 
lines relative to the continuum in composite spectra ( [Vanden Berk et al. 2001 ) . We use only H/3 and [Oni] AA4959, 5007 for 



the estimates and assume that no lines overlap except the BL and NL H/3. The source is at z s = 2. 

In this analysis we have ignored the noise associated with subtracting the sky and the absolute flux calibration. This 
might interfere with the measurement of the total flux from the NL, a, but should not seriously interfere with measurements 
of the ratios between NL and BL or between the continuum and the NL. We also ignored the fitting of the ratios between NLs 
which we do not think would affect our estimates by much. The line shape estimates are relatively statistically independent 
from the line ratio estimates so we also do not expect that having more complicated line shapes would greatly change our 
estimates. 



